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1. Introduction 

Consider a splierically symmetric spacetime filled by a massless scalar field The 
field \1/ evolves according to 

(7"^V„V,^' = , (1.1) 
where the geometry is determined by Einstein's field equations in the form 

Rab = ^'aa^'^f,* , (1.2) 

where = SttG. At first sight this problem would appear to be totally predictable and 
boring. Appearances are however deceptive. 

The system of equations (|1.1|) and ( p..2|) has been studied numerically, as a Cauchy 



problem in [0 and as a characteristic initial value problem in . In both cases the initial 
field can be specified in an arbitrary manner. Choptuik examined the following problem 
in [|1|. Suppose the initial data belong to a one-parameter set. Suppose further that for 
small values of the parameter p, the field is weak. One expects the evolution to be similar 



to the solution of (|1.1| ) on Minkowski spacetime; the field evaporates leaving an empty 
spacetime. Suppose that for large values of the parameter p, the initial fields are strong. 
Assuming some form of cosmic censorship hypothesis, we might expect a black hole to 
form around a spacetime singularity. By continuity there should be a critical parameter 
Pcrit separating the two types of evolution. 

Choptuik discovered numerically that nearly critical solutions appear to exhibit dis- 
cretely self-similar behavior — structure appears on ever finer scales [Q. (This was not 
discovered in previous numerical work because an adaptive grid algorithm is needed to 
resolve it.) If this behavior is real then the field and curvature can grow arbitrarily large 
close to the axis. Further, Choptuik exhibited in [|1| a scaling law. For supercritical evolu- 
tions [p > Pcrit) a black hole forms and he obtained a numerical law estimating the black 
hole mass 

Mbh = Ci{p - PcritV . (1.3) 

Here q depends on the nature of the initial data but 7 ~ 0.37 appears to be independent 
of the initial data. These results were confirmed in [0 , and similar results for other forms 
of matter have been established in . 

Discrete self-similarity is hard to treat analytically (see however Q), and so the- 
oreticians have looked for other forms of matter for which continuous self-similarity is 
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possible. Assuming continuous self-similarity for the critical solution reduces the problem 
to an ordinary differential equation. A perfect radiation fluid coupled to gravity was stud- 
ied in [Q, where it was numerically found that the critical solution is indeed continuously 
self-similar. Analytical studies of continuously self-similar solutions where the scaling sym- 
metry mixes with field symmetries have been made for the complex scalar field and for 
the axion/dilaton field 0, but it was not demonstrated that the continuously self-similar 
solution is indeed the critical solution. 

This paper reports on a numerical study of the axion/dilaton problem using the nu- 
merical techniques of f^. In §2 we describe the axion/dilaton problem in detail and §3 
describes the numerical procedure. Section 4 reviews the conjectured continuously self- 
similar critical solution. The main result of the paper is given in §5; there is a continu- 
ously self-similar solution which appears to be the attr actor at criticality. In §6 we consider 
slightly supercritical solutions and confirm both analytically and numerically the scaling 
law (|1.3|) . We estimate 

Tanalytic = 0.2641066 , 7numeric = 0.264 , (1.4) 

giving excellent agreement between theory and experiment. This critical exponent is how- 
ever very different from the value 7 ~ 0.37 found for a massless scalar field [Q, black body 
radiation 0] , and gravitational radiation . 



2. The axion/dilaton system 

We use the notation of 0. The matter used in this paper occurs in the 3-1-1- 
dimensional low-energy effective action of string theory. It consists of gravity coupled to a 
dilaton (j) and an axion p. The axion arises as the dual of the string three-form H, via 

Habc = -e'"t'eabcdd''p. (2.1) 

K 

The dilaton and axion can be combined into a single complex field 

r = 2p + ze"2'^. (2.2) 

(We have chosen to rescale p by a factor of 2 from the usual definition.) 
(presumably) irrelevant gauge fields, the effective action is 

' = 2^y^"^r"2a^i^J ' 
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Omitting the 



(2.3) 



where R is the scalar curvature. 

At the classical level, the model ( |2.3|) has an extra 5'L(2,1R) symmetry that involves 
only r and leaves the metric invariant. The symmetry acts on r as 

- - ^ • (2.4) 

where a,b,c,d & IR and ad — be = 1. This 5'L(2,IR) symmetry is presumably broken by 
quantum effects down to an SL{2, 2) symmetry, which is the popular S'-duality conjectured 
to be a non-perturbative symmetry of string theory. 
The equations of motions arising from (|2.3|) are 
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Rab - 77FirrT2 (^"'^^^'^ + daTdbT) = , (2.5a) 



V"V„r+— ^ = 0. (2.56) 

imr 

We will use two approaches to solve these equations. First, we will integrate them numer- 
ically given specified initial conditions. This will give us a range of numerical solutions, 
some ending with black holes, some ending with fiat space. Second, we will follow 
and assume that the exactly critical solution is continuously self-similar. The spacetime 
can then be derived analytically. We then compare the answers, and find extremely good 
agreement. 



3. The Numerical Approach 

The numerical integration is based on the approach described in . Throughout this 
paper, we assume spherical symmetry. As in [0, the coordinate system most suitable to 
the problem at hand is a double-null coordinate system. We take the spacetime metric to 
be 

ds"^ = —a^{u, v) du dv + r'^{u, v) dVL^. (3.1) 

Double-null coordinates are not unique. The remaining gauge freedom consists of redefining 
u and V hy u ^ f{u) and v g{v)-, where / and g are monotonic functions. We fix one 
of the functions by requiring u = v on the axis where r{u,v) = 0. The remaining gauge 
freedom will be fixed by the initial conditions. See 0] for further details. 
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With these coordinates, the equations of motion (|2.5| ) become 

r + rur^ + ^a^ = , (3.2a) 

auv/a - ttua^jc? + Vuvlr + (\)u(\)v + e^'^ PuPv = , (3.26) 

run - 2a^r^/a + r {(f>J + e^'^pj) = , (3.2c) 

r„„ - 2a,rja + r {(PJ + e^^ p^^) = , (3.2c^) 

r (puv + ru4>v + ry(j)u - 2r e^'^ puPv = , (3.2e) 
r Puv + T^uPv + f^Pu + 2r {(f) 

uPv + <?> vPu 

) = . (3.2/) 



It should be noted that when p = 0, the equations (|3.2[) reduce to a real scalar field coupled 



to gravity. This allows us to check much of the code. 

Since we shall be comparing results in different coordinate systems, we need to know 
what invariant quantities can be built from the above fields and metric. One quantity is 
the proper time T on the u — v axis, given by 

a(u, u) du . (3.3) 

All curvature invariants can be constructed from the following three building blocks. First, 
there is the Ricci scalar curvature 

^= (</'.<^. + e^W.) , (3.4) 

where we have used the equations of motion ( p.2|) . Two other curvature invariants exist, 
defined by 

^ V 32 ^^''^^''^ ~ " ~ '^"'^"^ ' ^^-^^ 



= J l^CabcdC'^''"' - = ^ (a' + 4rn^-) = 2m{u, v)/r^ , (3.6) 
V 12 D r^a^ 

where Cabcd is the Weyl tensor, and m(M, v) is the Hawking mass. It is simple to show, 

e.g., using boundary conditions given later in this section, that on axis, u = v, R2 = 0, 

and Rs = —R/6. 

We can convert the system of equations ( |3.2|) to a first order system by defining new 
variables 



^ = 4>y , ^ = (l)u , ^ = Pv , ^ = Pu , h = , h = Tu , e = a^/a 
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(3.7) 



In these new variables, fB.^) becomes 



r hu + hh + -a^ 

r h„ + hh + -a^ 
- 4 
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to which we adjoin 



/i„ -2eh + r [i)^ + e^-^^^j 



o^, — ae = 
r„ — /i = 
P. - ^ = 

4>y -tl; = 
















(3.8a) 

(3.86) 

(3.8c) 
(3.8d) 
(3.8e) 
(3.8/) 
(3.8^7) 
(3.8/i) 



(3.8z) 
(3.8i) 
(3.8A;) 
(3.80 



which foUow from eq. ( p.7| ). 

Next we write down the boundary conditions which apply on axis, u = v. Obviously, 
r = there, which implies + Vy = 0, i.e., h = —h. Eq. { 'i.Sc ) will be regular on axis if 
and only if a = 2/i (assuming a > 0). Now ( |3.8.q| ) implies ^ = ^ and (|3.8,f|) implies i/j = if; 
on axis. From the definitions (|3.7|) , these imply = <^r = on axis. It follows from a 
power series expansion (in powers of r) that = on axis. 

Suppose that the fields e, i/;, and ^ are known as functions of f on a surface u = uq = 
constant. We can integrate (|3.8/|) as an ordinary differential equation in v for (p. The initial 
data come from the boundary condition (pr = on axis. Eq. ( |3.8/c|) for p and ( |3.8i| ) for 
a are treated in identical fashion. Next, eqs. ( |3.8d| ), (|3.87"|) are integrated as simultaneous 
equations for r and h, with r = 0, /i = a/2on axis. Eq. ( p.86|) with data h = —h on axis 
gives /i as a function of v. Finally we integrate ( |3.8/] ), ( |3.8/j|) with data ip = ip, C = 
on axis, to obtain tp and ^ as functions of v. Thus we now have all quantities on the 
surface u = uq. Next we determine 6,-0, and ^ on the surface u = uq + Au = constant by 
integrating ( |3.8c| ), ( p.8^ ) and (|3.8g| ). Returning to the start of this paragraph, we continue 
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the integration process until either a singularity forms or the fields have evaporated away 
leaving fiat spacetime. 

In order to start the integration we need to give initial data e, '0, and ^ on the initial 
surface u = 0. The remaining coordinate freedom is equivalent to the choice of e (or a) on 
the surface u = 0. This is arbitrary. The initial fields and ^ determine the physics. 

The u-integration is done using an explicit finite difference algorithm. The v- 
integration is done using an implicit (midpoint) algorithm. However, if the equations are 
integrated in the manner described above, they are linear and so the method is effectively 
explicit. We operate exclusively at a Courant number of 1, and the implicit algorithm 
provides the necessary stability. There is no need for any artificial viscosity. 

During the evolution fine structure appears on ever-decreasing scales. We use the 
adaptive mesh algorithm of Berger and Oliger ||] as modified in 0. This enables us to 
insert extra grid points (again with a Courant number of 1) where and when necessary to 
resolve the fine detail. 

4. The analytic critical solution 

It was conjectured in that the precisely critical solution of the axion/dilaton system 
coupled to gravity is continuously self-similar. The continuously self-similar solution is 
described in detail there, but we review it here. Continuous self-similarity means that 
there exists a homothetic Killing vector ( that satisfies 

jCc9ab = '^9ab ■ (4.1) 

The homothetic symmetry mixes with the SL{2, IR) symmetry, changing r by an infinites- 
imal 5'L(2,1R) transformation 

Cc^T = ao + aiT + Q!2r^ , (4.2) 

with ai G M. 

The standard approach to spherically symmetric self-similar solutions is to assume 
the metric has the form* 

ds^ = (1 + u{t, r)) (-62(t, r) dt^ + dr^) + dQ^ . (4.3) 

* Ideally, one would use double-null coordinates as in the preceding section. For unknown 
technical reasons, we have been unable to repeat the calculations of in a it — f coordinate 
system. 
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The time coordinate t is chosen so that the singular point occurs at t = 0, and the metric 
is regular for t < 0. The homothetic vector is ^ = tdt + rdr- One then defines a new 
scale invariant coordinate z = —r/t. In terms of the z — t coordinates, the metric is 



ds^ = {l + u{z)) [{z^ -b^{z)) dt^ + 2ztdtdz + t^ dz^] + z^t^ dQ^ , (4.4) 

and C, becomes simply C, = tdt- After fixing unimportant constants, the axion/dilaton field 
has the form 

rit.z)=^ ]-tZ^^;} , (4.5) 



where u; is a real constant, and f{z) is a complex function. We will sometimes write 
f{z) = fmiz)e'^f''^^\ where fm and fa are real functions. 

Details of the construction of the solution can be found in 0. The r = axis 
corresponds to z = and is a singular point of the equations of motion. One can choose 
b{z = 0) = 1 for t < 0. The past light cone of the origin corresponds to the line where 
b{z+) = 2+, and is also a singular point of the equations of motion. Requiring regularity 
along the r = axis (for t < 0) and regularity at z^ uniquely fixes the solution in terms of 
w, z+, 1/(0)1, and \fiz+)\, with 

u = 1.17695272200 ± 0.00000000270 , (4.6a) 

z+ = 2.60909347510 ± 0.00000000216 , (4.66) 

1/(0)1 = 0.892555411872 ± 0.000000000224 , (4.6c) 

1/(2+) I = 0.364210875022 ± 0.000000000760 . (4.6d) 

We can now construct the same invariant quantities as in the u — v coordinates. On 
the r = axis, u{z = 0) = and 6(0) = 1, so the proper time is simply T = t. The 
curvature invariants are 



R 



t'^z{l + u{z)) 



2a;|/(z)P/:(z) _u{zl 



' 2t%{z){l + u{z)){l-\f{zWr' 

R, = }^ (4 9) 

^ t^z^il + uiz)) ■ ^ ' 
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We can use regularity on the r = axis to determine that 

2u;^|/(0)P 53.38 
^(^ = 0) = - r2(i_|;(o)|2)2 = ' (4-10) 

and -R2(^ = 0) = 0, -R3(^ = 0) = — i?/6, as expected. On the nuU cone b{z+) = z+, there is 
no proper time. Define A = (l + (u;^ — 2)\f{z+)\'^ + |/(2+)|^). Using regularity, we find 

^ t2^2A(l + (^2_2)|/(,^)|4 + |/(,J|8) ' ^^-^^^^ 

o , ^_ U^\f{z,)?{l-\f{z,)mi+\f{z.)?) 

4t2;,2A(l + (^2_2)|/(^J|4 + |/(^J|8) ' ^^-^^^^ 

^3(..) = ^^^1^ . (4.11c) 

To compare with the numerical solution, we will find it more convenient to consider just 
the ratios, or 

0.7866, (4.12a) 



Rs{z^) l + {u;^-2)\f{z^r + \f{z. 



R{z,) 2|/(^,)P(l-|/(zJP)2 
R2{z,) {l-\f{z^r)\l + \f{z^)?) 



Rs{z^) Ml + {u;^-2)\f{z^)\^ + \f{z^)\^) 



-4.958 , (4.126) 
= -0.1587 . (4.12c) 



5. Comparison of approaches at the critical solution 

As in P] , we can tune the initial conditions to be very close to the threshold of black 
hole formation. This should, in principle, result in a solution very close to the exactly 
critical solution of the previous section. This behavior is most clearly displayed on the 
r = axis, where we expect the numerical solution to resemble most closely the exactly 
critical solution. 

The scalar curvature on axis as a function of proper time is shown in fig. 1. The initial 
conditions chosen for the numerical solution shown contain much more dilaton than axion, 
but the late time behavior is independent of initial conditions. The x-axis is — log(Tci-it —T) 
(all logarithms in this paper have base 10), so the system is evolving from left to right. 
For the analytical solution, Tcrit = 0. For the numerical solutions, we are free to adjust 
Tcrit to most closely match the analytical solution. The y-axis is log(— i?). Three lines are 
shown in fig. 1. The solid line is a slightly sub-critical numerical evolution. The curvature 
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log(-i?) 




2 2.5 

- iog(rcHt - 

Fig. 1: The scalar curvature as a function of proper time on axis. The soHd hne is 
a sHghtly sub-critical numerical solution. The dashed line is a slightly super-critical 



numerical solution. The dotted line is the theoretical prediction (4.10). 



becomes strong on axis, but eventually dissipates (when T > Tcrit). The dashed line is a 
slightly super-critical solution. It becomes a black hole, as can be seen by its upturn for 
late T. The dotted line is the continuously self-similar solution ( [4.10| ). 

As expected, on the left side of fig. 1, the numerical and analytical solutions dis- 
agree. This is the region where the numerical solutions depends on the initial values of 
the fields. As the evolution progresses, the numerical solutions approach the analytical 
solution, becoming almost indistinguishable for a while. This is strong confirmation that 
the continuously self-similar solution is indeed the attractor. Eventually, the numerical 
solutions diverge from each other and from the analytical solution. 

We compare the solutions on the past lightcone of the singular point in fig. 2, fig. 3 
and fig. 4. From the analytic point of view, this is straightforward. Eq. (|4.12|) tells us 
that the ratios of the curvature invariants should approach constants. Numerically, it 
is much trickier. The past lightcone corresponds to a line of constant v = VcAt- One 
must choose fcrit carefully because any error will grow as it approaches the r = axis. 
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4 4.5 

— log r 

Fig. 2: The ratio R2/R as a function of — logr along a line of constant v. The solid 
line is a slightly subcritical numerical solution. The dotted line in the theoretical 
prediction ( [i.l2a|) . As r ^ 0, R2/R 0. 



Unfortunately, this is precisely where the numerical solution should be approaching the 
analytical solution. To further complicate matters, since the lightcone is a null curve, 
there is no invariant definition of distance along the curve. In the figures, we have plotted 
the curvature ratios as a function of — log r. The function r becomes multivalued for 
supercritical evolutions (which is why r is a bad choice as a coordinate), so we show a 
slightly subcritical numerical evolution in the figures. We see that the numerical solution 
(shown as a solid line) approaches the analytical line (the dotted line) , remains there for a 
while, and then diverges from the analytical solution. When r = 0, we expect i?2 = and 
i?3 = — i?/6, and these are indeed the values seen for very small r. The transition from the 



z = values (|4.12|) to the r = values is fairly abrupt. The transition is at finite r because 
of errors in choosing f crit and because the numerical solution is not precisely critical. Note, 
for example, that i23(logr = —4) ^ 10^ and is growing quickly near logr = —4. Thus the 
fact that the ratios are roughly constant and close to the analytical prediction ([4.12| ) is an 
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1.5 2 2.5 3 3.5 4 4.5 5 5.5 



— log r 

Fig. 3: The ratio Rs/R as a function of — logr along a line of constant v. The solid 
line is a slightly subcritical numerical solution. The dotted line in the theoretical 
prediction ( ^T^ . As r ^ 0, R3/R -R/6. 

excellent indication that the solution is indeed the continuously self-similar solution given 
in §4. 

6. Moving away from criticality 

One of the most interesting findings in [|l|] was the discovery of the scaling law ( |1.3| ) 
for the mass of black holes away from criticality. The striking result was that 7 was 
completely independent of initial conditions, and was thus dubbed the "universal" critical 
exponent. For the case of a real scalar field coupled to gravity, numerical work has given 
7 = 0.374 [00]. Previous numerical work with other types of matter gives similar values for 

7. A study of axisymmetric gravitational wave collapse 0] indicates 7 ~ 0.36. Spherically 
symmetric perfect radiation fiuid Q gives 7 ~ 0.36. This coincidence of numbers led 
people to suspect that 7 might be a universal exponent, independent of the specific type 
of matter coupled to gravity. 
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3.5 
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log r 



4.5 



5.5 



Fig. 4: The ratio R2 / R3 as a function of — log r along a line of constant v. The solid 
line is a slightly subcritical numerical solution. The dotted line in the theoretical 
prediction ( p^ . As r ^ 0, R2/ R3 ^ 0. 



The universality of 7 has been drawn into question by analytical arguments. Starting 
with a known critical solution, one can use renormalization group techniques to calculate 7. 
This was first carried out for the perfect radiation fiuid in |]^, which found 7 = 0.3558, in 
good agreement with [Q, but only fairly close to the value for a real scalar field. Because 
of the numerical uncertainties in |^ and uncertainties in the accuracy of the analytic 
arguments, it is unclear whether these results indicate conclusively that 7 depends on the 
matter content. A more general class of perfect fiuids with p = kp was studied analytically 
in |T^ under the assumption that the critical solution is continuously self-similar. It 



was found that 7 depends strongly on k. However, no numerical work has been done for 
/c 7^ 1/3, and it is not clear that the continuously self-similar solution really is the attractor 
at criticality. 



The case of a complex scalar field has been studied in [|Tl|. They also perturbed 
away from a continuously self-similar solution, and found 7 = 0.3871, which seems to be 
distinctly different from the perfect radiation fiuid value. However, they also found that 
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the continuously self-similar solution is unstable, and conjecture that it is an attractor of 
codimension three. The discretely self-similar solution (the same critical solution as for 
the real scalar field case) has codimension one, and therefore dominates. Thus numerically 
it should be difficult to find the continuously self-similar solution. We have numerically 
evolved the complex scalar field with a variety of initial conditions and in every case have 
foimd critical behavior of the discretely self-similar type. Thus, no conclusive evidence has 
so far been presented that 7 is not universal. 

Given an exactly critical analytical solution, one can perturb the critical solution to 
find the critical exponent 7 as follows ||^ . Let h be any function of the analytical solution, 
such as 6 or /. Perturb away from the critical solution 

ht{z, t) = hssiz) + e|tr"/lpert(-2) , (6.1) 

where /iss(^) is the critical solution, e is a small number, k is a constant, and /ipert(^) 
depends only on z. Replacing ht{z, t) in the equations of motion and keeping terms to first 
order in e gives an eigenvalue equation for k. This eigenvalue equation can in principle 
have a number of possible solutions for k. The solution with the largest value of Re k will 
cause the fastest growing perturbation in (|6.1j), and is called the "most relevant mode." 
The critical exponent is given by P, p!Oi|TT| 

7 = • (6.2) 

Re K 

We have carried out a perturbation analysis of the axion/dilaton critical solution 
described in §4. In principle, k can have an imaginary component, but this greatly com- 
plicates the analysis, and we found it to be unnecessary in this instance. The result for 
perturbation of the continuously self-similar solution is 

Tanalytic = 0.2641066 . (6.3) 

This value is quite different from the critical exponents found in previous numerical studies, 
and therefore should be easy to distinguish numerically. Our analysis does not rule out 
other possible values of k, but we have seen no numerical evidence for other unstable 
modes. 
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log(p-pcrit) 

Fig. 5: Mass as a function of p — pait for the axion/dilaton system. 



We show in fig. 5 the numerical results for a large number of supercritical evolutions 
along a one-parameter family of initial conditions parameterized by p. The mass is defined 
as in 0. Pick a constant vq and determine the point [uq^vq) where the apparent horizon 
intersects vq. The mass of the black hole is defined as Mbh = rn{uo, vq), where m{u,v) is 
the Hawking mass defined in ( p.6| ). As can be seen from fig. 5, log Mbh is a linear function 
of \og{p — Pcrit) until one is far away from the critical solution. Omitting the final seven 
points in fig. 5 where the curve is no longer linear, we obtain a slope for fig. 5 of 

Tnumeric = 0.264 , (6.4) 

which agrees nearly perfectly with the theoretical prediction ( |6.3| ). As in the case of the 
real scalar field, 7numeric seems to be independent of the initial data, as long as the axion 
in not identically zero. This is further strong evidence that the critical solution is indeed 
continuously self-similar, and that the perturbation expansion gives the right answer. 

7. Conclusions and unanswered questions 
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Fig. 6: Points along the critical line in the axion/dilaton system. The a;-axis 
measures the strength of the axion initial data, and the y-axis measures the strength 
of the dilaton initial data. 



We depict a two-dimensional subspace of the infinite-dimensional space of initial con- 
ditions in fig. 6. The x-axis is the (normalized) amplitude of the initial axion field, and 
the y-axis is the (normalized) amplitude of the initial dilaton field. The points on the 
graph represent critical solutions. Not surprisingly, the critical surface has codimension 
one. When the axion field is initially identically zero, it remains zero, and the problem 
reduces to a real scalar field coupled to gravity. The critical value for that solution is 
the point on the y-axis, and we know from that the critical solution there is discretely 
self-similar. The conjecture of is that the points on the critical line with nonzero axion 
are continuously self-similar. This raises the question: is the discretely self-similar solution 
simply an artifact of special initial conditions? In other words, is the discretely self-similar 
solution an isolated point on the critical line, or is there an open set on the critical line 
near the y-axis where the critical solution is discretely self-similar? This is a difficult ques- 
tion to answer numerically, because near the y-axis, the axion is quite small and is thus 
liable to be swamped by numerical noise. The solution in fig. 1 represents a point quite 
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near the y-axis, and is clearly continuously self- similar. By moving along the critical line 
even closer to the y-axis, we push the region in fig. 1 that depends on initial conditions 
farther to the right. Eventually, this will meet the point where the sub- and super-critical 
solutions diverge from the analytical solution. Thus we cannot expect to move arbitrarily 
close to the y-axis in a numerical evolution. Except for points at which we believe the 
numerical accuracy is breaking down, every point near the y-axis that we have checked is 
continuously self-similar. This implies that the discretely self-similar solution is unstable, 
but analytical work is necessary to confirm this. 

The axion/dilaton matter system described in this paper is physically well motivated. 
However, if the discretely self-similar solution is simply the result of non-generic initial 
conditions, might the continuously self-similar solution also be the result of overly restricted 
matter? For instance, string theory contains numerous gauge fields in the low-energy 
effective action which we have omitted in writing (|2.3|) . Including these fields might result 
in yet another type of critical solution. If we hope to use critical solutions to probe the 
small-scale structure of spacetime, we should in principle use the most general solution 
allowed in our favorite theory of quantum gravity. 
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